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The existence of pentaquarks, namely baryonic states made up of four quarks and one antiquark, 
became questionable, because the candidates, i.e. the O"*" peak, are seen in certain reactions, i.e. p+p 
collisions, but not in others, i.e. e'^'e" annihilations. In this paper, we estimate the production of 
©"'"(1540) and H (1860) in e"'"e~ annihilations at different energies using Fermi statistical model as 
originally proposed in its microcanonical form. The results is compared with that from pp collisions 
. . . at SPS and RHIC energies. We find that, if pentaquark states exist, the production is highly possible 

\0 \ in e+e" annihilations. For example, at LEP energy ^=91.2 GeV, both 6+ (1540) and H (1860) 

, yield more than in pp collisions at SPS and RHIC energy. 
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I. INTRODUCTION 



Pentaquark is the name for baryons made up of four quarks and one antiquark. Recent advances in theoretical and 
experimental work led to the observation of pentaquark candidates by some experiments, however, other experiments 
, reported no observation from their search. For more details, see the review of the experimental evidence on pentaquarks 
and critical discussion The question of the existence of pentaquarks is raised again, due to the non-observation 
in so many experiments. To answer the question, it is important to assess the cross sections in different processes. 

The candidate 0+ peak is seen in different reactions namely of 7+ A, 11+ A, p+p, K+Xe, e+d, e+p, K+Xe0|3l- 
All of those reactions involved at least one baryon in the initial state. But other experiments, e.g. e+e": BalDar, 
^ Belle, Bes, LEP experiments 0, pp: CDF, DO pA:E690, jp: FOCUS, pA: HERA-B, ep: Zeus (for the 0°) pi+ ^LiD: 
^ COMPASS, Hadronic Z decays: LEP, tt, K, p on A: HyperCP, 77: L3, 7r,p,I] on p: SELEX, pA: SPHINX, E'A: 
Qh, WA89 K^p: LASS 5], did not observe those candidates. Because the experiments such as e"*"e~ did not observe those 
• candidates, it is generally believed a non-zero initial baryon number is essential to the pentaquark production, i.e. 
D ' pp collisions or collisions involved with nuclei, as we can see the collision types of the positive reports. To check this, 
we calculate in this paper the pentaquark production in e"'"e~ annihilations. A theoretical comparison of pentaquark 
^ [ production between p-l-p0 and e+e" processes may be helpful to the question if pentaquarks exist. 

Why we choose the microcanonical approach to calculate? Traditionally, the hadron production in e+e" scattering 
at high energy is treated as a two-stage process. First e+e~ — > -y*/Z qq is calculated using perturbative theory. 
Then the qq system produces hadrons phenomenologically based on string fragmentation or cluster fragmentation. 
Pentaquark states are exotic, hard to treated within the frame of conventional string models. In this situation 
statistical approaches may be of great help0,H]- It was Hagedorn who introduced statistical methods into the strong 
interaction physics in order to calculate the momentum spectra of the produced particles and the production of 
strange particles^]. Later, after statistical models have been successfully applied to relativistic heavy ion collisions 
[To|. Becattini and Heinz came back to the statistical description of elementary reactions. 
According to the situations of obeying conservation laws, statistical models are classified into four ensembles: 

* microcanonical: both, material conservation laws (Q, B, S, C, • • • ) and motional conservation laws (E, ~j) , J , 
• • •), hold exactly. 

* canonical: material conservation laws hold exactly, but motional conservation laws hold on average (a temperature 
is introduced). 

* grand-canonical: both material conservation laws and motional conservation laws hold on the average (temperature 
and chemical potentials introduced). 

* mixed ensemble: for example, to estimate particle production in heavy ion collisions, so-called partially canonical 
and partially grand-canonical models are employed where strangeness conserves strictly, but the conservation of net 
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charge hold on the average. 

What one expects is that the microcanonical ensemble must be used for very small volumes, i.e. the systems created 
by the collisions between elementary particles. For intermediate volumes the canonical ensemble should be a good 
approximation, while for very large volumes, i.e. the systems created by heavy ion collisions, the grand-canonical 
ensemble can be employed. Therefore we take a microcanonical approach to calculate the production of pentaquark 
states from electron-positron annihilations. 



II. THE APPROACH 

The idea was originally proposed by Fermi in its microcanonical form0 and realized with Markov chain 
technique|3|- We calculate the hadron production in e''"e~ annihilations at a given energy as a statistical decay 
of a cluster which carries net quark contents Q = (iV„ — N^, Nd — N^, Ng — Ng) = (0, 0, 0). The cluster is charactered 
by three parameters, cluster energy (mass) E, volume V and strangeness suppression factor 7^. We assume that 
hadron production from the cluster is dominated by the n-body phase space. More precisely, the probability of the 
cluster hadronization into a configuration K — {hi,pi; . . . ; hn,Pn} of hadrons hi with four momenta pi is given by 
the microcanonical partition function il{K) of an ideal, relativistic gas of the n hadrons 12], 

"^""^ = (2^ n^^^" n — l[d'P^^iE-^s.)Sm)SQ^^,^, (1) 

^ ' i=l aes i=l 



with Ei — ^ nif + p^ being the energy, and pi the 3-momentum of particle i. 

The term Sq^j^q- ensures flavour conservation; is the flavour vector of hadron i. The symbol S represents the set 
of hadron species considered: we take S to contain the pseudoscalar and vector mesons {n, K,ri,r]' , p, K* ,uj, cj)) and 
the lowest spin-^ and spin-| baryons {N, A, E, S, A, S*, S*, f2) and the corresponding antibaryons. ria is the number 
of hadrons of species a, and gi is the degeneracy of particle i. 

It's well known that strangeness will be overpopulated if the hadron production is purely determined by the 71-body 
phase space. The common treatment is to introduce so-called strangeness suppression factor js with < 7s < 1. So 
we employ this factor as well in Eq. (1), and the index is the number of (anti)strange components in the final-state 
particle i, i.e. for Kaons, A, E and their antibaryons, Si — 1; for 0-meson, 2 and S, si = 2; for f2 and 17, si = 3. 

Similar to the previous work[^, we add the pentaquark states 9"''(1540) and S(1860) into S. The O"*" has quark 
contents [uudds]. The 2(1860) can be S~~ [ddssw], S^[dssuu] , 'EP[ussd(]\ or S+[MMssd]. The spin of pentaquark 
states can not be determined by experiments yet, and it is generally accepted pentaquark states are spin-i particles, 
so we take a degeneracy factor g = 2. As for strangeness suppression, = 1 for 8+ and Si = 2 for S(1860). 

The high dimensional phase space integral is verified via constructing Markov chains of hadron configurations K. 
The Metropolis algorithm provides random configurations K according to the weight, the corresponding microcanon- 
ical partition function fl(K). All possible random configurations are included. 

Working with Markov chains one has to worry about two kinds of convergences: the number of iterations per chain 
must be sufficiently big (essentially a multiple of the so-called auto-correlation time), otherwise the method is simply 
wrong. Secondly, the number of simulated chains must be sufficient large, to obtain the desired statistical accuracy. 
Questions related to the auto-correlation time have been studied in detail in earlier publications 0, so that the error 
due to auto-correlations can be neglected. The statistical error will be treated carefully at the result section. 

In addition to checking auto-correlations and statistical errors, also physics cross checks have been performed: a 
comparison of our Monte Carlo method with canonical method gives a good agreement when the systems have 
big volumes, i.e. 50fm'^ and high energies, i.e. 10 times the mass of observed particles. 

We generate randomly configurations K according to the probability distribution fl[K). The Monte Carlo technique 
allows to calculate mean values of an observable as 

A = Y,AiK)niK)/Y,niK'), (2) 

K K' 

where ^ means summation over all possible configurations and integration over the pi variables. A{K) is some 
observable assigned to each configuration, as for example the An multiplicity Mh{K) of hadrons of species h present 
in K. Since Mh depends on E, V and 7s, we usually write Mh{E, V,js)- 

The cluster parameters energy E, volume V and strange factor 7s are chosen so that the cluster decay reproduces 
best the measured multiplicity of the selected hadrons in e+e~ annihilations at a given energy of ^/s. This is achieved 
by minimizing x^' 

X-(i.,l^,7.)^E'''^"^^'^^^"^'^''''^'"^^^' (3) 
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TABLE L Cluster parameters for e+e annihilations at diflterent energies determined by the yields of tx^ ,p, , K and their 
antiparticles. 





£(GcV) 


V(fm=') 


Is 


xVdof 


=29 GeV 


10.8±0.8 


62±24 


0.64±0.11 


1.97/5 


^/i =35 GeV 


12.4±0.8 


70±20 


0.71±0.14 


4.97/5 


y/s =91.2 GeV 


17.6±0.4 


102±12 


0.61±0.07 


2.92/5 



where Moxp,j('\/s)i ^^^d Uj are the experimentally measured multiplicity and its error of some selected hadron species 
j in e^e~ collisions at an energy of ^/s. 

With the cluster parameters determined by the selected hadrons, we can predict the multiplicity of any particles 
included in the hadron species list, i.e. pentaquark states from e+e~ collisions at the energy of -^s. 

III. THE RESULTS 

We calculate pentaquark production from e+e" collisions at a/s = 29,35,91.2 GeV. This is realized by two steps: 
first, we have to determine the cluster parameters energy i?, volume V and strange factor 7^ for each collision energy 
•v/i; then, we use the obtained cluster parameters to calculate pentaquark production. 

For the first step: with data taken from we select j = 7r+,p, if+,A and their antiparticles to determine the 
cluster parameters energy E, volume V and strange factor 73. The available experimental yields of other particles 
are used to check the reliability of our results. For each set of parameters {E, V, js), the multiplicities of selected 
hadrons are calculated from 10,000 random hadron configurations so that the relative statistical errors of the selected 
hadrons are within 0.03. Then the best-chosen set of parameters is obtained by minimizing ^-s in Eq.(3). Fig^ 
displays the results of our fit procedure in comparison with the experimental data at different energies. Here the error 
bars of theoretical yields are statistical. 

Fig.l tells that the best-chosen cluster parameters with the experimental multiplicities of tt^ ,p, , A from e+e^ 
collisions at the energy of ^/s = 29 GeV can also reproduce very well the yields of other particles such as 7r°, p'^, t], cj), S~ . 
The less well reproduced particles are r]' , E*"*" and 0, spin-1 meson and spin-| baryons. 

At the energy of ^/s = 35 GeV, the best-chosen cluster parameters with the experimental yields of tt^ ,p, , A 
can also reproduce very well the yields of other particles such as tt^ ^ , rj, K*^ ,E.~ . 

At the energy of y/s = 91.2 GeV, the best-chosen cluster parameters with the experimental yields of 7r+,p, A 
from e+e^ collisions can well reproduce the yields of other particles such as vr*^, p^, 77, cj, if*, E, S^. The less well 
reproduced particles are 77', A"*""*" , E* , S* and fl, spin-1 meson and spin-| baryons. 

The large list of reproduced particle proves that the cluster parameters, energy E, volume V and strange factor 
7s , determined by 4 particle yields of 7r+ , p, , A can be reliably used to estimate the yields of particles, especially 
for spin-0 mesons and spin-i baryons. As it is accepted that the spin of pentaquark states is our estimation on 
pentaquark yields does not suffer from the problem with spin-3/2 baryon yields. 

In Tab. Qlthe best-chosen set of parameters (i?, V 7s), for e''"e~ collisions at the energy of ^/s =29, 35 and 91.2 GeV 
are collected. The small x^/dof indicates a good fit quality. 

We find that with the increase of e+e~ energy y^, the created cluster has a higher energy E and a bigger volume 
V . The strange factor 7s does not change much with the increase of e+e" energy ^/s. 

One may feel hard to understand why E = y/s does not hold, instead, the created cluster has an energy E much 
smaller than the e"''e~ energy y^. Let's imagine an expanding fireball (mainly longitudinal/thrust) created in e+e~ 
annihilation. In this case, there is a lot of collective kinetic energy. The parameter E means the sum of the energies of 
volume elements in their proper frames, in other words their invariant masses. So this effective mass is much smaller 
than the mass of the total system, which is of course e+e" energy y^. If we only consider total yields, we do not 
need to specify the details of the collective expansion. The price is that this simple 3-parameter model cannot make 
any statements about transverse momentum spectra or rapidity spectra. 

The volume we find here is much bigger than that in Fermi's model 0, where the volume is estimated according 
to the colliding particles with the Lorentz contact. In our case, the volume V is the size of the collision system at 
the moment of hadronization. After e+e~ annihilates, the collision system gets many quark pair production and 
expanding with time. So the size is much much bigger than that of an electron and a positron. 

The strangeness suppression factor 7s changes very little with the collision energy ^/s. The value is around 0.7, bigger 
than what we found in pp collisions. This is also consistent with the result from canonical fitting by Becattini[ll|. 

With the cluster parameters obtained above, we can calculate the production of all hadrons in the hadron list. 
Pentaquark states have been included in the hadron list, so we can get the 47r yields of pentaquark states from e+e^ 
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FIG. 1: (Color online) A comparison between hadron production from cluster decay and from e^e~ annihilation experiments 
at 29 GeV, 35 GeV, and 91.2 GeV. 



annihilations at the energy of 29, 35, and 91.2 GeV. Since the yields of pentaquark states are much lower than ordinary 
hadrons, we calculate their yields from the average of 1,000,000 random configurations with Eq.(2). In Fig. |21we plot 
the result as 47r yields of pentaquark states versus the center mass system energy y^, and compare with that from 
pp collisions at SPS energy Eiab=158 GeV and RHIC energy Ecms=200 GeV. The error bar is a combination of the 
statistical Monte Carlo evaluation error with the influence of the flt error of the parameters. The two errors are at 
the same magnitude, and the statistical one is smaller. 

The production of pentaquark states, i.e. and S , from e+e^ annihilations increase with the collision energy. 
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FIG. 2: A comparison of the 4tv yields of O"*" (left) and H (right) from e"'"e~ annihilations and from pp collisions at0. 

In pp collisions, we can see the similar behavior for all pentaquark states except 0+. The production of 8+ in pp 
collisions is relatively high at a few GeV, decreases slightly, then increases at teens GeV and keep increasing for higher 
collision energiesj^. The special behavior of 0+ in pp collisions at low energies is due to the proton excitation, which 
does not apply in e+e" collisions. 

The yields of 8+ from e+e" collisions at the energies we study are at the same magnitude of that from pp collisions. 
The 8+ production from e+e" at LEP energy ^/s =91.2 GeV is even higher than the yields from pp collisions at both 
SPS and RHIC energy. 

For S , we get obviously higher production from e"'"e~ collisions than from pp collisions. This can be understood: 

1. The average mass of clusters created in e^c^ collisions is much bigger than that in pp collisions, i.e. the clusters 
have average mass 17.6 GeV for e+e" at y/s =91.2 GeV, and mass 7.3 GeV for pp collisions at SPS and mass 
16.15 GeV at RHIC energy. The cluster mass is most sensitive parameter for particle production. Certainly cluster 
with big mass can produce more particles. 

2. The strangeness suppression factor is bigger in e'^e^ collisions. S contains two strange quarks, and gets squared 
strangeness suppression, which is about 0.1 in pp collisions and about 0.4 in e+e~ collisions. 

The particle ratios 0'^ /p and S~~(1860)/S~ are of interest and have been discussed by several different approaches 
in heavy collisions and pp collisions. Here we study the ratios in e^e~ collisions at different energies. 

We get the /p about 0.03 in e+e" annihilations, while it is much smaller, about 0.007, in pp collisions with 
the microcanonical approach '6^. We recognize the big difference can be caused by the st rang eness suppression factor, 
0.6 ^ 0.7 in e+e" annihilations, but 0.33 in pp collisions. The grandcanonical ensemble [l5j gives about 0.06 for this 
ratio in heavy ion collisions, which is even higher. This occurs when the strange chemical potential /i^ ~ 0, which 
corresponds to 7^ = 1, an even bigger strangeness suppression factor. 

We get higher S (1860)/S~ ratio, about 0.1, in e+e~ annihilations, which is 0.02 in pp collisions with the 
microcanonical approach 6], and 0.01 in heavy ion collisions obtained by the grandcanonical results'l65. 

IV. DISCUSSION AND CONCLUSION 

We estimate the pentaquark production, i.e. 8+ and S , from e+e~ annihilations at the energy of y/s—29, 35 
and 91.2 GeV using Fermi statistical model as originally proposed in its microcanonical form. We obtain increasing 
production of pentaquark states with the increase of collision energy. Comparing with the previous work for pp 
collisions, we find that the yields of 8"*", from e+e" at the above mentioned energies are at the same magnitude as 
those from pp collisions at SPS and RHIC energies. The yields of 8+ from e+e" at LEP energy y/s—91.2 GeV is 
higher than the yields from pp collisions at SPS and RHIC energies. 

From our estimation, the production of S from e"*'e~ collisions is obviously higher than that from pp collisions at 
SPS and RHIC energies. That's very different from the experimental report - the observation of S was reported for 
the first time in the SPS experiment, by NA49 collaboration's!; but none of the four collaborations in LEP experiment 
get the observation, though they did search. 

Theoretically we conclude that Fermi statistical model pentaquark production is quite high from e+e^ annihilations 
if pentaquark states do exist due to a very high energy clusters are created when the e+e" annihilate. The energy 
parameter E=17.6 GeV for e+e" at y/s =91.2 GeV, while E=16.2 GeV for pp collisions at y/s =200 GeV. The 
average proper mass of clusters created in pp collisions at RHIC energy is lower, indicates a lot of energy is taken 
away collectively, due to the existence of the leading particles. Initial baryons are not necessary for pentaquark 
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production. The clusters with big masses created from any kind of high energy colhsions can provide a rather high 
yield. 
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